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We give a simple and efficient process for generating a quantum superposition of states which 
form a discrete approximation of any efficiently integrable (such as log concave) probability density 
functions. 
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It is well known that probabilistic algorithms are some- 
times more powerful (and often more practical!) than 
their deterministic counterparts in classical computing. 
Such algorithms generically make use of some source of 
random bits. It has been shown Q however, that giving 
a quantum computer access to random bits does not, in 
fact, further empower it. 

There has been some attention recently to the problem 
of using quantum computers to generate certain proba- 
bility distributions, in particular by looking at the prop- 
erties of quantum random walks . 

In this short note, we will look at the question of 
whether, given a certain probability distribution {p{\, 
we can efficiently create a quantum superposition of the 
form 



\^({pi}))=J2VPi\i), 



(1) 



where the \i) are an orthonormal set of states. We are, of 
course, particularly interested in the case where the index 
i ranges over some exponentially large set of possibilities 
N. 

We have not solved the general problem of when a 
state of the form (I) can be efficiently created. However 
we will show here that if an efficient classical algorithm 
exists to integrate a certain probability distribution p(x), 
then we can efficiently create a state of the form (I) 
for the discretized version {p{\ of p(x). A well known 
set of probability density functions which are efficiently 
integrable by monte carlo methods are log- concave dis- 
tributions. (A log-concave distribution is one for which 

9 gx^^ < 0-) Most important statistical distributions 
(including the exponential and normal families) are log- 
concave. The task of sampling and/or integrating such 
distributions frequently arises for problems of determin- 
ing the (multidimensional) volumes of solids || . 

Clearly if an efficient classical algorithm exists, then we 
can certainly create a quantum computer in the mixed 
state p = ^2iPi\i){i\- Our interest in creating a state of 
the form (1) is thus not simply to be able to measure the 
state |z) with probability pt. Rather we hope that gen- 
crating the distribution in a coherent fashion will allow 



further (uniquely quantum) processing of this distribu- 
tion. At the end we discuss some possibilities along these 
lines. 

We will illustrate the procedure considering only a 
probability distribution over a single random variable, 
although the method can be easily extended to multi- 
variable distributions. 

We let n — logN, where N is the total number of 
points over which we wish to discretize this distribution. 
Imagine the distribution is divided into some number 2™ 
of regions, and that we already have a m qubit state 
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where, as usual, the integer i, is interpreted in binary as 
the tensor product sequence of qubit states. Here p^ is 
the probability for the random variable x to lie region i; 
thus Pg™^ is the probability for x to lie in the far leftmost 
region, p^™ 1 ' is the probability for it to lie in the region 
adjacent to this and so on. Our goal is to show that we 
can now efficiently subdivide these 2 m regions to yield a 
2™ l+1 region discretization of p{x). That is we wish to 
add one qubit to the state (2), such that we achieve the 
evolution 
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(3) 



where a (resp. (3) is the probability for x to lie in the left 
(resp. right) half of region i. If we can achieve such an 
evolution, then the new state of our m + 1 qubits is 



\lpm+l) = Yl 



(m+1) 



(4) 
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This process is then repeated until m = n, i.e. until 
we have created the desired superposition (1) over all 
N = 2™ states. If such an evolution is possible, then 
clearly it is efficient. 

To show how an evolution of the form (3) can be 
achieved, we first define x l L and x l R to be the left and 
right boundaries of region i. We then define the function 
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J x i L 2 p(x)dx 
J*i R p{x)dx 



2 



f(i) is simply the probability that, given x lies in region 
i, it also lies in the left half of this region. Since J p(x)dx 
is efficiently computable classically, we can take an an- 
cilla register initially in the state |0 . . . 0), and construct 
a circuit which efficiently performs the computation 

^\i)\0...0)^^\i)\9 i ), (5) 

where 6i = arccos \J f{i). We now perform a controlled 
rotation of angle 6i on the m + l'th qubit: 

v5F|i>|«>-> ^\i)\e i )(co S 9 i \0} + sin^|l)) (6) 

and we then uncompute the register containing \9i) to 
leave us in a state of the form (3) as desired. Note that 
the efficient classical algorithm to perform the integra- 
tion may well be probabilistic, however it can always be 
implemented on a quantum computer in the required co- 
herent manner by using an ancilla register initially pre- 
pared containing the string of random bits to be used in 
the computation. 

Having seen that we can create a state of the form (1), 
we now give some examples of why creating such a state 
might be interesting. 

(i) Non-uniform priors for quantum searching The 
quantum search algorithm [|4| is typically envisioned as a 
search over a large number of possible solutions to a prob- 
lem, where every potential solution is a-priori equally 
likely. In many instances however, we might believe that 
the correct solution is more likely to be found in a cer- 
tain region of state space than another. Thus we might 
wish to input a non-uniform prior distribution into the 
search algorithm. If this prior distribution is log-concave 
(e.g. gaussian, poissonian etc) then the procedure out- 
lined above can be used to generate it. 



(ii) Sampling of non-log- concave distributions A state 
of the form (1) is clearly of little value if we were to just 
measure the qubits and obtain outcome \i) with proba- 
bility pi, since such sampling is possible classically. The 
fundamental difference between a quantum superposition 
and a classical sample however, is quantum interference. 
In particular, we may perform unitary transformations 
on this state to take it to some other state of the form 

3 

Each of the new probabilities qj can in some sense be 
functions of the full set of Pi's, and we may obtain "in- 
terference" between the probabilities that has no classi- 
cal analogue. For example, a simple Walsh-Hadamard 
transformation on each of the the qubits produces the 
new distribution 

* = ^(£(-ifVft) 2 (8) 

i 

and this distribution is not log-concave. 

(iii) Estimating the magnitude of a fourier component. 
A state of the form (1) can be fourier transformed ef- 
ficiently ||. If we are interested in the magnitude of a 
partciular fourier component, then techniques exist for 
the process of amplitude estimation, and these proceed 
quadratically faster than their classical counterparts. 



Acknowledgments 

This research was supported by the NSA & ARO under 
contract No. DAAG55-98-C-0040. 



[1] D. Aharonov, A. Kitaev and N. Nisan. Proceedings of 30th pg. 156. 

STOC, 20 (1998). [4] L. Grover, Phys. Rev. Lett. 79,325 (1997). 

[2] D. Aharonov, A. Ambainis, J. Kempe and U. Vazirani, [5] P. Shor, SIAM J. Comp., 26, 1484 (1997). 

Proceedings of 33th STOC, 20 (2001). 
[3] D. Applegate and R. Kannan, Proceedings of STOC '91, 



